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Abstract. We present the results of quantitative spectral analysis of sub- 
dwarf O stars (sdO) obtained from the SDSS database. By visual inspection, 
112 sdO stars could be identified in a set of spectra of photometrically selected 
faint blue stars. Fitting the Balmer, He I and He II lines to state-of-the-art non- 
LTE model spectra, we derived their effective temperatures, surface gravities 
and helium abundances. We find the helium-enriched sdO stars to rally in a 
smah intervaU of Teft = 40 000 K ... 50 000 K and log g = 5.5 . . . 6.0, whereas the 
helium-deficient ones are showing a wide variance both in temperature and log f/. 
A puzzling feature is a significant number of helium rich subdwarfs below the he- 
lium main sequence. These stars would not be able to burn helium in their cores 
and pose a serious problem to subdwarf evolution theories. We conclude, that 
the sdB stars are linked with the helium-deficient sdO stars, i.e. the sdBs are 
predecessors for the helium-deficient sdOs. To explain the helium-enriched sdO 
stars in their entirety, we regard two scenarios as most promising: The merging 
of two helium-core white dwarfs and the late hot fiashers. To investigate this 
further, we started an NLTE analysis of carbon abundances. So far, we find 
carbon slightly enriched above solar, but no trends can be seen yet. 



1. Introduction 

We distinguish spectroscopically between the cooler sdB stars (Tgff < 40 000 K) 
showing He I li nes and the hot sd O stars (Tefj > 40 000 K), with Hell lines. 
Recent work of IStroer et alj ( 20071 ) suggests a strong correlation of the hehum 



abundance with carbon/nitrogen line strengths. A physically more meaningful 
classification into 'helium-enriched sdO' for stars with supersolar helium abun- 
dance and 'helium-deficient sdO' for subsolar abundances is suggested. 

While the sdBs form a homogenous group, the sdOs show a wide spread 
in temperature, gravity and helium abundance. Only a few stars are helium- 
deficient but the helium-enriched sdO stars are represented in great number. 
Canonical theories see an evolutionary link between sdO and sdB stars, though 
the mechanism driving the hydrogen-rich sdB atmosphere into a helium-rich sdO 
in the course of stellar evolution is hard to explain. 

We are interested in the sdO stars, particulary their evolutionary status, 
their origins and possible evolutionary links between them and the coo l er sdB 
stars. Recent work in this field has been carried out by IStroer et al.l ( 20071 ) 



using the high resolution spectra of 46 sdO stars from the ESO SPY survey 
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( Napiwotzki et al. 200ll ). We will take advantage of the huge database of the 



ongoing Sloan Digital Sky Survey which provides spectra of many different ob- 
jects, among them a number of sdOs. By comparing and combining our results 
with the results from iStroer et akl (2003), the sample size will be large enough 
to test rivalling theories for their origin and evolution 



2. Quantitative Analysis of SDSS Spectra 

From the database of the 5th data release (DR5) of SDSS, we selected all point 
sources within the colour box [u — g) < 0.4 and {g — r) < 0.1. We classified more 
than 8 000 spectra by visual inspection, mostly white dwarfs, about 500 sdB 
stars and 112 sdOs. After sorting out spectra with spectral features indicative 
of a cool companion star or of too low signal-to-noise, we carried out a spectral 
analysis of the remaining 86 of our programme stars. A model atmosphere fit 
was conducted in order to obtain the effective temperature, surface gravity and 
helium abundance. 

2.1. Model Atmosphere Fitting 
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Figure 1. Examples of NLTE model fits to sdO spectra from SDSS. Left 
hand side: a helium-deficient sdO. Right hand side: a helium-enriched sdO. 
The observed spectra are displayed as histograms. Lines that are not fitted 
but only displayed for inconsistency checks are plotted as dashed lines. 



To determine Teg, log 51 and log(^^^) of the stars, we fitted synthetic model 
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spect ra to the observed ones, using a x -minimization procedure ([Napiwotzki 



19991 ) to derive all three parameters as well as the radial velocity simultane- 



ously. The synthetic non -LTE models were constructed using the PR 02 code 
(IWerner Dreizleij[l999l ') with a new temperature correction technique ([Dreizlerl 



2OO3I ) and include partially line blanketing. It includes H and He atoms only. 



The model grid ranges in temperature from 30 000 K to 100 000 K and from 
4.8 dex to 6.4 dex for log^ with log(^^-) = —4 . . . -|- 3. No extrapolation beyond 
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the model grid was allowed. Due to the inconsistencies repo rted for the Hn,-lin e 



of some sdB stars, which might be caused by a stellar wind ( Heber et al.ll2003l ). 
this line was never fitted to the models, but it was always plotted to search for 
peculiarities. to and He II 4859 A and 4339 A were always fitted, but we 
found He II 4100 A often too weak to provide a reasonable fit. The ionization 
equilibrium of He II and He I is very sensible to temperature, therefore we always 
included He 1 4472 A, even in its absence. 

2.2. Considering the Errors 

The fit program provided us with mere statistical errors, they range from AT^s ~ 
150... 3 000 K, depending on the signal-to-noise ratio. Fortunately, six stars 
were observed twice (see Table. [1]), which we can use to estimate the quality of 
our fits. The differences in parameters obtained from spectrum A and spectrum 
B lie within the given statistical errors, which we therefore conclude to be robust. 



Table 1. Two exposures had been taken from six of our programme stars. 
We use this opportunity to check the reliability of the errors. 



SDSS designation Toff log.g log(^) 

[K] [cms-2] 



SDSSJ 


00:51:07.01+00:42:32.5 


38234+330 


5 


73+0 


06 


-1.07 






38207+305 


5 


72+0 


05 


-0.94 


SDSSJ 


11:14:38.57-00:40:24.1 


57714+3268 


5 


56+0 


19 


-0.94 






55730+2292 


5 


56+0 


15 


-0.86 


SDSSJ 


16:37:02.78-01:13:51.7 


45547+338 


5 


77+0 


08 


+3 






45860+353 


5 


70+0 


08 


+3 


SDSSJ 


17:00:45.67+60:43:08.4 


47404+1870 


5 


92+0 


16 


+2.26 






49537+2911 


6 


09+0 


17 


+2.17 


SDSSJ 


23:35:41.47+00:02:19.4 


69906+2234 


5 


45+0 


08 


+ 1.10 






72870+1341 


5 


44+0 


11 


+ 1.18 


SDSSJ 


23:39:13.99+13:42:14.2 


48592+1540 


5 


82+0 


15 


+0.92 






48914+3233 


5 


66+0 


18 


+ 1.19 



2.3. Results and Comparison to the SPY Sample 

We now compare our SDSS sample to the 46 stars of the SPY sample. We find 
that the ratio of helium-enriched stars to helium-deficient ones is higher in SDSS 
than in SPY (68:18 vs. 33:13). In addition, the SDSS sample contains more very 
hot stars (Tefj > 60 000 K) than the SPY sample. Both findings are probably 
due to selection effects, e.g. the limiting magnitude of the surveys, which is 
considerably fainter for SDSS than for SPY. Besides that, the distributions of 
the SDSS and the SPY sample in the Tefj-log g plane are very similar, as can be 
seen from Fig. [5J We will therefore merge the samples for further discussion. 

None of the programme stars is located on the EHB, meaning that either 
they are post-EHB stars, most likely evolved from the sdB stars, or they do 
not have any link to the EHB at all and are just evolving through this re- 
gion in the HR diagram. The helium-enriched sdOs are spread over Tgfj = 



4 



Hirsch, Heher, and O'Toole 



... *■ . 




70 60 50 

T,„ / 1000 



/ 



/ 



60 50 
T.„ / 1000 



Figure 2. Left: Tcff-log(-^) diagram. The solid horizontal hne represents 
the solar helium abundance. Right: Tcs-logg diagram. Also given is the 
helium main sequence (HeZAMS) as well as the EHB, defined by the location 
at zero age and terminal age. Filled symbols represent helium-enriched stars, 
while the open symbols are helium-deficient ones. Squares are stars from 
SDSS, triangles are stars from SPY. 



40 000 K . . . 74 000 K and log (7 = 5.2 . . . 6.3, but a concentration of them within 
the small temperature interval T^q = 40 000K . . . 50 000K becomes apparent, 
only a few stars are found to lie at hotter temperatures. The hydrogen-rich 
sdOs on the other hand shun this region. 

Note that both the SDSS and the SPY sample contains stars lying below 
the helium main sequence, all but one are helium-enriched. These stars cannot 
sustain stable helium burning in their cores and therefore it is difficult to explain 
their origin. 



3. Stellar Evolution Tests Using the Combined Samples 

Two different aproaches to hot subdwarf star evolution are presented in the 
literature: (i) binary star evolution and (ii) single star evolution. Both are trying 
to explain the high mass loss the stars have to suffer from before reaching the 
EHB. In the context of binary star scenarios, this is achieved by mass exchange 
in sufficiently close systems, whereas single star scenarios have to invoke physical 
mechanisms within the star itself, like unusually strong stellar wind, mixing, etc. 



3.1. Binary Star Evolution 

Han et aTl ( 20021 . 20031 ) did an extensive binary population synthesis and found 



three formation channels relevant for the sdB/sdO population: stable Roche 
lobe overflow, common envelope ejection and a merging of two helium-core white 
dwarfs. While the first two scenarios result in relatively cool stars preferentially 
lying between 30 000 K and 40 000 K, our sdO stars are very unlikely to have 
followed the latter formation channel as they are mostly hotter. 
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For sdO stars, especially the helium-enriched ones, we regard the forma- 
tion via a merging of two helium WDs as most promising. Short period binary 
WDs will lose orbital energy through gravitational waves. With shrinking sep- 
aration, the less massive object will even tually be disrupted an d accreted onto 
its companion, leading to helium ignition. ISaio k Jeffervl tOOCf) argue, that this 



merger product will result in a helium buring subdwarf showing an atmosphere 
enriched in CNO-processed matter. This scenario therefore can explain these 
extremely helium- enriched sdOs showing s trong nitrogen lines in their atmo- 
spheres. However, iGourgouliatos Sz Jeffervl (|2006) find that, under the assump- 
tion of total angular momentum conservation, He+He WD merger do rotate 
faster than breakup velocity. A mechanism that enables the star to get rid of 
its angular momentum still has to be found. 

However, there is another possible origin for hot subdwarf stars. The sdB 
HD 188112 is reported to have a mass of 0.23 Mq, base d on a Hipparcos p arallax 
measurement, too low to sustain core helium burning ( Heber et al. 20031 ). Such 
a star may have formed if mass transfer occured on the first red giant branch 
(RGB) and before the onset of core helium bu rning. 

Apropriate tracks have been calculated bv lDriebe et al.l ( 19981 ). who evolved 
a 1 Mq main sequence star as one component of a binary system through the 
RGB stage without ignition of helium to a helium star and eventually to a 
helium core white dwarf, under the assumption of mass transfer to the compan- 
ion. The evolution of such stars with masses between 0.414 Mq and 0.300 Mq 
would lead straight through the observed sdO population, down to below the 
HeZAMS. However, how compelling these finding s may be, they suffer from a 
fundamental discrepancy: iNapiwotzki et ah ( 20041 ) find that the fraction of ra- 
dial velocity variable stars among the helium-enriched sdOs is 4 % at most. Also, 
the predicted helium abundance is much lower than for the helium-enriched sdOs 
suitable for this scenario and within the relevant mass range, no enrichment of 
the atmosphere with processed matter from the interior is predicted. 



3.2. Single Star Evolution 

The canonical evolution theory considers sdO stars to be the progeny of the 
sdB stars irrespective of their helium content. Thus a star reaching the EHB 
will evolve from the EHB towards the sdO regime with higher temperatures and 
eventually to the white dwarf coolin g curve. Such tracks have been calculated by 
Dorman. Rood Sz O'Connell ( 19931 ). They started with a star at the tip of the 



RGB and calculated the minimum core mass needed to ignite helium burning, 
then added a small envelope. The needed core mass varies slightly and ranges 
from 0.464 Mq for supersolar to 0.495 Mq for subsolar metallicity, while the 
envelope mass is very small (Menv < 0.002 Mq). Still, the question for the cause 
of the mass loss remains unanswered. In Fig. [3] we compare three post-EHB 
evolution tracks to the observed distribution. They follow the evolution from 
the zero age horizontal branch, where core helium burning starts, to the terminal 
age horizontal branch, and beyond, for stars with 0.471, 0.473 and 0.475 solar 
mass cores. Because of the small envelopes, the post EHB evolution hydrogen 
shell burning is insignificant and therefore the star does not climb the AGB 
(ABG manque stars), but evolves through higher temperatures until eventually 
it settles on the white dwarf cooling curve. These post-EHB tracks cover the sdB 
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Figure 3. Left: Comparison of our sample to three post-EHB tracks 
(iDorman. Rood fc O'Connell [T993h for 0.475 A/q, 0.473 M^t, and 0-471 M,t, 
(top to bottom). In the upper left two post AGB tracks (iSchonberneilll983l ) 
for 0.565 Mq and 0.546 Mq (left to right) are shown as well. They may ex- 
plain a few of the very hot sdO stars. Right: A track of a late hot flasher after 
it has settled onto the EHB is plotted (|Sweigart.l997. : .Mochlcr et al.,,2002i) . 
Symbols as in Fig. ^ 



regime as well as the lower gravity helium-deficient sdO stars at Tgfj > 40 000 K. 
This is consistent with the helium-deficient sdO stars being the descendents of 
the sdB stars. 

However, some of the hottest stars and in particular the helium-enriched 
stars which cluster at Tgfj = 45 000 K do not fit into this scheme. The very hot 
sdO stars may be explained as post-AGB star s on their way to th e white dwarf 
cooling curve. Two such tracks calculated bv ISchonberne 2 (|1983I ) are shown in 
Fig. [3l But due to the very short evolutionary timescales (~ 30 000 yrs from 
ejection of the planetary nebula to pre-white dwarf settling) for the post-AGB 
stars, we expect only a few such stars in our sample. 

For the helium-enriched sdOs on the other hand, we have to find another 
mechanism that explains both their distribution as well as the peculiar car- 
bon/nitr o gen l ine s trengths in connect ion with the high helium abundances. 
SweigartI (|l997l ) andU rown et al.l (|200lh analyzed a scenario for stars called late 



hot flasher. They argue, that fast rotation of RGB stars could lead to mixing of 
helium into the envelope, resulting in a higher RGB peak luminosity and there- 
fore to higher mass loss rates on the RGB. This will delay the helium flash until 
the star has left the RGB and is already descending on the white dwarf cool- 
ing curve. The delayed flash will induce mixing which will transport hydrogen 
into the helium burning core where it is burnt, resulting in a helium-enriched 
star with Teff ~ 40 000 K on or near the helium main sequence and enriched in 
carbon. 

We have plotted an evolutionary track resulting from a late hot flasher in 
Fig. [3l Also this track fails to reproduce the observed distribution, though the 
very late hot flasher scenario can explain stars below the helium main sequence. 
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4. Carbon Abundances from High-Resolution Spectra (SPY) 

It is clear now, that the hehum abundance and the T^s-logg diagram alone 
cannot answer the question as to the sdOs' origin. As different scenarios predict 
different abundance patterns in particular for carbon and nitrogen, we started 
with a carbon abundance analysis of high resolution spectra of the sdO stars 
from SPY. The new models include C III and C IV in addition to hydrogen and 
helium. They were calculated and fitted to the observed spectrum using the 
same methods as descibed above. 

4.1. Preliminary Results 

So far we have investigated seven helium-enriched sdOs. All stars are slightly 
enriched in carbon, with log( ./^^ ) = —2.6... — 3.0, compared to the solar 

value of log( j^*^ ^ ) = —3.5. However, when fitting the carbon lines, two stars 
showed a significant line broadening indicating projected rotational velocities of 
25kms~^ (see Fig. HD and 30kms~^, which is surprising as most sdB stars do 
not show 't;rotsini > lOkms"^ (see Geier et al., these proceedings). Though 
low number statistics, this might be considered another hint pointing towards 
the assumption that sdB stars and helium-enriched sdO stars are two different 
populations that formed by different evolutionary paths. 




° 3600 3610 4180 4190 4150 4160 5800 5810 

wavelength in A 



Figure 4. Line fits to C III 3609.0, 3609.6 A, C III 4186.9 A, C III 4152.4, 
4156.5, 4162.8A and CIV 5801.3, 5812.0A (from left to right). Top panels: 
No rotational broadening, bottom panel: Wrotsini = 25kms^^. 



5. Summary and Conclusion 

We have fitted synthetic NLTE model spectra to 87 sdO spectra from the SDSS 
database, using a technique. While the helium-enriched stars cluster in a 
relative small regime around Tefj = 45 000 K, the helium-deficient sdO stars are 
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widely spread. Worth noting is a significant number of lielium-deficient sdOs 
below the helium main sequence, a fact that cannot be explained by canonical 
evolution. 

Preliminary results of a carbon abundance analysis shows slightly supersolar 
values for seven helium-enriched sdOs. The higher Vrot sin i we found in two of 
them compared to the very slowly rotating sdBs might be an additional support 
for the merger scenario. 

By comparing the stars' distribution in the T^s-logg diagram to several 
evolutionary paths, we conclude that most likely the helium-deficient sdO stars 
are the progenies of the sdB stars. The helium-enriched sdO stars still pose 
a problem: The late hot flasher scenario as well as the white dwarf mergers 
qualify as formation channels as they predict an enrichment of helium and 
carbon/nitrogen. The merger scenario is attractive because it can also natu- 
rally explain the absence of radial velocity variables (close binaries) amongst 
the helium-enriched sdOs. 
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